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Abstract 


This paper presents the performance analysis of a planar solid-oxide fuel cell (SOFC) under direct internal reforming conditions. A detailed 
solid-oxide fuel cell model is used to study the influences of various operating parameters on cell performance. Significant differences in efficiency 
and power density are observed for isothermal and adiabatic operational regimes. The influence of air number, specific catalyst area, anode 
thickness, steam to carbon (s/c) ratio of the inlet fuel, and extend of pre-reforming on cell performance is analyzed. In all cases except for the case 
of pre-reformed fuel, adiabatic operation results in lower performance compared to isothermal operation. It is further discussed that, though direct 
internal reforming may lead to cost reduction and increased efficiency by effective utilization of waste heat, the efficiency of the fuel cell itself is 
higher for pre-reformed fuel compared to non-reformed fuel. Furthermore, criteria for the choice of optimal operating conditions for cell stacks 


operating under direct internal reforming conditions are discussed. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Given a cell geometry, the performance of a SOFC strongly 
depends on the operating conditions and the inlet fuel composi- 
tion. Since SOFCs offer a wide range of operating possibilities, 
identifying the effect of operating conditions such as air flow 
rate and inlet fuel composition on the efficiency and power den- 
sity are critical for optimal operation of SOFCs. For example, 
extent of pre-reforming will affect fuel utilization and power 
density. On efficiency maps, operating in the region of high effi- 
ciency may result in a very low power density requiring large 
cell volume to deliver the required power output, while operat- 
ing at high power density can lead to low efficiency [1]. Hence, 
the SOFC often has to be operated in a region of reasonable 
efficiency and power density. On the other hand, material cost 
can be optimized by judicious choice of catalyst loading and cell 
geometry. For example, an anode supported cell is a preferable 
choice for hydrocarbon fuels. However, the amount of catalyst 
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loading and anode thickness for optimal performance are design 
parameters [2]. 

In an earlier study, Zhu and Kee presented a maximum 
efficiency analysis using a thermodynamic model [1]. They con- 
cluded that the maximum possible efficiency is independent 
of the characteristics of membrane electrode assembly and the 
internal polarization losses. A number of publications report 
exergy analysis of SOFC systems [3-5]. Unlike these models, the 
study presented here is based on a very detailed SOFC model and 
focuses on a unit cell. In a previous study, we presented a numer- 
ical model to simulate a planar cell under isothermal conditions 
[6]. The study incorporated detailed models for heterogeneous 
chemistry using an elementary step reaction mechanism for the 
steam reforming of CH4 on nickel catalysts. Furthermore, the 
electrochemistry was modeled using a modified Butler-Volmer 
setting based on single-step electron transfer reaction. In the 
present study, the model presented in [6] is extended to non- 
isothermal conditions in order to analyze the effects of the 
operating variables on cell performance. Though it may be pos- 
sible to run a cell under isothermal conditions in a laboratory 
environment, the scenario is completely different for a techni- 
cal system. In general, each channel in a SOFC stack behaves 
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Nomenclature 

Ac area of cross section of flow channels (m?) 

As specific area (m7!) 

Bg permeability (m7) 

Cy specific heat (Jkg~! K7!) 

dp particle diameter (m) 

D diffusivity (m? s7!) 

Dh hydraulic diameter (m) 

Ecell cell voltage (V) 

Eh thermoneutral voltage (V) 

F Faraday number (C mol!) 

Gz Graetz number 

h heat transfer coefficient (J m~* K7! s~!); specific 
enthalpy (Jkg~!) 

He channel height (m) 

AH enthalpy of formation (Jmol-') 

i current density (A cm~?) 

Jk species flux (mol m~? s7!) 

k thermal conductivity (J m`! s7! K7!) 

Kg number of gas-phase species 

m mass flux (kg s- 1) 

n number of charge transferred 

n molar flow rate (mol s7!) 

Nu Nusselt number 

Omin minimum O; required (mol m) 

p pressure (Pa) 

Pe MEA perimeter (m) 

Pr Prandtl number 

Q Source term (J (s m3)~!) 

R gas constant (J mol! K7!) 

Re Reynolds number 

5 molar production rate (mol m~? s~!, 
mol m~? s7!) 

AS entropy change (J K7!) 

t time (s) 

T temperature (K) 

v velocity (m sh 

W average molecular weight (kg mol~!) 

We molecular weight of kth species (kg mol!) 

[X] concentration (mol m”) 

Y mass fraction 

Greek symbols 

r surface site density (mol cm~?) 

ô Kronecker symbol 

€ emissivity 

n efficiency 

Nu fuel utilization 

0 surface coverage 

À air number 

u viscosity (kg m7! s7!) 

p density (kg m~?) 

oO co-ordination number; conductivity (S cm!) 

T tortuosity 

{0) porosity 


Subscripts 

c channel 

e electrolyte/electrode 
f fluid 

I interconnect 

k species index 


differently. For example, the temperature profile developing at 
the outer most channels in the stack can be quite different 
from the one located at the center of the stack. Nevertheless, 
it is rational to assume that the extremes of the stack behav- 
ior can be well understood by isothermal and adiabatic cases. 
Therefore, calculations are performed for both isothermal and 
adiabatic conditions. However, running detailed models like the 
one reported in this paper on stack level would require enormous 
computational time. Therefore, a single cell is analyzed in the 
present study. 


2. Model description 


In a previous publication, we reported an elliptical model to 
describe the transport and chemical processes within the cell 
[7]. However, in this paper we discuss a parabolic model to 
describe the operation of a planar SOFC, which reduces the com- 
putational effort of the mass and heat transport calculation. The 
following sections encompass the modeling framework adopted 
in this work. A schematic representation of the model geometry 
considered in this study is shown in Fig. 1. 


2.1. Mass and heat transport 


2.1.1. Channel flow 

Flow through fuel and air channels can be assumed to be 
one-dimensional and laminar. Then the plug flow equation for 
species continuity 


aloYr) _ _ ƏloruYr) Pr Pe 


JkWk, k=1,..., Kg, (1) 


ot əz Ac 


__--—— Interconnect 


_Anode 
Pr 


~#——- Electrolyte 


Cathode 


Fig. 1. A schematic representation of planar cell under consideration. 
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and the mass balance equation 


dr ahu) 
ot dz 


Kg 
+ u) Pode We, (2) 
k=1 


can be applied. Here pr is the fluid density, Y; the species mass 
fraction of species k, u the velocity, Pẹ the perimeter associ- 
ated with the membrane electrode assembly, m the mass flux, z 
the axial position, ¢ the time, and W% is the species molecular 
weight. Assuming constant pressure in the channels, the density 
is calculated from the ideal gas equation 


pW = prRT. (3) 


In Eqs. (1) and (2), Jg is the flux at the electrode channel 
interface and is calculated using the dusty gas model (DGM) 
(Eq. (11)). The temperature 7; in the flow channels is determined 
from the heat balance equation 


KprC ptT) KuprCpT) h h 
= T; — T. Ty — Tẹ). 
at 3z H., f e) + aC I f) 
(4) 


Here Cpr is the specific heat capacity of the fluid; Tf, Te, and Ty 
are the temperature of fluid stream, electrode, and interconnect, 
respectively. The first term on the right-hand side of Eq. (4) 
represents the transport of energy due to the bulk fluid flow, 
second term represents the heat transfer from the channel to the 
electrode structure, and third term represents the heat transfer 
from the interconnect into the flow channels. The heat transfer 
coefficient h is evaluated from the Nusselt number 


_ ADp 
=a 
where Dp is the hydraulic diameter and k is the thermal conduc- 


tivity of the fluid. The Nusselt number is expressed empirically 
as [8]: 


1000\ 70-5386 6.7275 
Nu = 3.095 + 8.933| —— exp | — (6) 
Gz Gz 


where the Graetz number, Gz, is given in terms of the Reynolds 
number, Re, and Prandtl number, Pr, as 


Nu (5) 


Dh 

Gz = — Re Pr, (7) 
Z 

where z is the axial position along the channel. Re and Pr are 


defined as 


(8) 


2.1.2. Porous media transport 

Species transport through the porous media is assumed to be 
one-dimensional along the thickness of the porous structure and 
given by 


0 Y; a( Jg W; 
(dprYk) u (Jk Wx) 4 WA. (9) 
ot oy 


Here, & is the porosity, 5, the heterogeneous molar production 
rate of the chemical species k, y the independent spatial variable 


along the thickness, and A, is the specific catalyst area available 
for surface reactions. Total density of the fluid within the porous 
structure can be calculated from 


K Ke 
adr) <A Si We) SQ, 

= Wz As. 10 

F 3 ay + Dik kAs (10) 


In the above equations, the fluxes J; are evaluated using the 
dusty gas model (DGM). According to DGM, the net species 
molar flux is given by [6,9-11] 


Ke Ke 


PPM ix | B 
Jk =— DPEMVIX a Ze ay 
k > kl [Xi] + 3 Di, nn” 


(1) 


The first term on the right-hand side of Eq. (11) represents the 
diffusive flux and the second term represents the viscous flux. 
a are defined as DGM diffusion coefficients given as [6] 
DDM — AO}, (12) 


where the elements of H matrix are 


1 X; Xk 
= rn a (13) 
k,kn j#k kj kl 

The permeability B, in Eq. (11) is given by the 

Kozeny-Carman relationship [12]: 
eds 

~ 7211 p) 
Here dp is the particle diameter and t is the tortuosity. Assuming 


the reaction heat is released on the solid surface the heat balance 
equation within the porous structure can be written as 


g (14) 


Kg 
= V (keV T) 5 WrAshk (15) 
k=1 


(pC pT) 
ot 


Here, p is the density of the solid structure, Cp the specific heat 
capacity of the solid structure, and kerr is the effective thermal 
conductivity. The first term on the right-hand side represents 
heat transfer due to conduction and the second term repre- 
sents the heat release due to heterogeneous reactions within 
the electrodes. Three source terms, Qr, On, and Qe, enter as 
the boundary condition for Eq. (15). The radiative heat transfer 
between the interconnect and the outer most discretized cell in 
the porous electrode is given by 


de 1 o (T} - T*) a6) 
"~~ dy | (/er) + (/e) — 1]’ 


where ez and e represent the emissivity of the interconnect and 
the electrode, respectively. The convective heat transfer between 
the fluid in the channel and the outer most discretized cell in the 
porous electrode is given by 


h 
Qn = yT =T). (17) 
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dy appearing in Eqs. (16) and (17) results from the finite volume 
integration over the discretized cells and must have a finite value. 
Qe arises from the electrochemical heat release at the three- 
phase boundary (Eq. (19)). The global charge transfer reaction 
at the three-phase interface can be written as 


H + 0?- > HO+2e”, (18) 


however, only a part of the energy change during the reaction is 
released as heat, that amounts to 


pa, (19) 


Assuming the electrolyte to be one-dimensional the heat bal- 
ance equation is reduced to 


d(eCpT) af aT\ 7? 
=— |k ; 20 
ot əz ( T) = Oe (20) 


The last term in the above equation represents the ohmic heat- 
ing due to ion transport in the electrolyte. Though the electrodes 
and interconnects also contribute to ohmic heating due to elec- 
tronic resistance, their contribution is minor compared to the 
ionic resistance offered by the electrolytes and is neglected in 
the present study. 


2.1.3. Interconnect 
The heat balance equation for the interconnect can be sum- 
marized as 


(oC _ 9 / ƏT h 
ot ~ az (x əz ) T m” T) 
4 nd 
a ee (21) 
Hi | (1/€&) + A/e) - 1 


Here, Cpr is the specific heat capacity, and p is the density of 
the interconnect material. The first term on the right-hand side 
represents the conduction within the interconnect, second and 
third terms represent the heat transferred to/from the channel to 
interconnect and the radiation transfer to/from the interconnect 
to the electrode structure. 

Eqs. (15), (20), and (21) require boundary conditions at z = 0 
and z = L given by 
oT oT 


= 0, = 0. (22) 


dz z=0 az z=L 


Modeling radiative heat transfer in SOFCs is a complex process 
and the radiative transport within the semitransparent electrodes 
and electrolyte as well as surface to surface heat transfer must 
be accounted for an accurate calculation. However, an exact 
knowledge of phenomenological properties like absorption coef- 
ficient, refractive index, scattering coefficient, emissivity, and 
reflectivity is not available and hence represents an obstacle in 
modeling radiative heat transport in SOFC electrodes and elec- 
trolyte. A detailed discussion of radiation heat transfer in SOFC 
is given elsewhere [13-16]. In an analysis of spectral radiation 
in SOFC electrodes, Damm and Fedorov [14] proved that the 
radiation effects in SOFC electrodes are minimal and can safely 


be neglected. Therefore, in the present work, we only consider 
surface to surface radiation. Since the planar geometry consid- 
ered here is of high aspect ratio (L/d > 50) the surfaces can be 
treated as black bodies with unity emissivity [13]. 


2.2. Heterogeneous chemistry and electrochemistry 


The anode chemistry model adopted here excerpts the poten- 
tial of elementary heterogeneous mechanism for the steam 
reforming of methane on Ni catalysts. The modeling frame- 
work of heterogeneous chemistry is described in our previous 
reports [6,7,17]. The relation between surface coverages 6, and 
the surface concentrations [X+] is given by 


_ [Xlo 
= 


and the temporal variations of surface coverages are given by 


OK (23) 


dé SkOk 
d Tr 
where 5; is the surface reaction rate, og the co-ordination number 
(number of sites required for a species for adsorption), I” the 
surface site density, and K, is the number of surface species 
[7,17,6]: 

Assuming H? as the only electrochemically active species and 
neglecting the contribution by CO, the electrochemistry model 
is implemented using a modified Butler-Volmer setting based 
on single-step electron transfer mechanism, and is described 
in detail in our previous report [6]. Neglecting the contribu- 
tion by CO is a safe assumption due to the low electrochemical 
activity of CO compared to H when both are present together 
in the system [18]. The activation energy and pre-exponential 
factor in the expression for the exchange current density for 
hydrogen oxidation and oxygen reduction are taken from [7]. 
Since the chemistry and electrochemistry models are discussed 
thoroughly in our previous reports, they are not a subject of 
further discussion here. However, we ascertain here that the het- 
erogeneous chemistry mechanism is not validated for carbon 
deposition conditions. 


, k=Kgt+1,...,Kg+ Ks, (24) 


2.3. Efficiency model 


There are many discussions in the literature on the efficiency 
analysis of fuel cells [1,19-21]. The efficiency of a cell is usu- 
ally expressed in terms of cell voltage, fuel utilization, and the 
thermoneutral voltage [22] 


Eeell 
Ein 


N= Mu (25) 
Here, nu is the fuel utilization, Ece the operating cell voltage, 
and Ep is the thermoneutral voltage. The fuel utilization, 


_ Wp idx 


: (26) 
fuen F 


Nu 
depends on the operating voltage, fuel and air flow rates, avail- 
able specific area for internal reforming, and the cell geometry 
[22]. In the above equation, W is the width associated with the 
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Table 1 
Cell geometry and operating conditions for various cases under consideration 
Case Length (cm) Fuel inlet condition 800 °C Air number Comments 
CH4% H20% u (cms™!) 
(a) 10 40 60 30 Variable — 
(b) 10 40 60 30 1 Variable anode thickness 
(c) 10 - - 30 3 Variable s/c ratio 
(d) 10 40 60 30 1 - 
(e) 5 - - 30 5 Air inlet at 600° C and various levels of pre-reforming 


MEA and /1fye1 is the molar flow rate of the fuel. The thermoneu- 
tral voltage is defined as 


AH 
nF’ 
where AH is the enthalpy of formation, and n is the total number 
of electrons transferred. 


En = (27) 


3. Computational procedure 


In order to solve the equation systems they are first cast into 
a finite volume form. The channels, interconnects, and elec- 
trolyte are treated as one-dimensional, and the reactor geometry 
is discretized in the axial direction into 200 cells. The anode is 
discretized into 25 cells along the thickness and the cathode into 
10 cells. Eqs. (1)-(4), (9), (10), (15), (20), (21), and (24) form a 
system of coupled non-linear equations which can be written in 
the residual form as 


F(®) = 0, (28) 
where the vector ® is given by 
® = [(Tic, (Y, m, T), O, p, T, Mar,» + 


(T)e, (Y, P; Tot .. 


+g (Y, P, T, an 


<, Œ, p, Deo Dei] - 
(29) 


Here the indexes (ic),(f),(a,1),(a,n),(e),(c,1),(c,m),(a) stand 
for interconnect, fuel channel, first discretized cell in the anode, 
nth discretized cell in the anode, electrolyte, first discretized cell 
in the cathode, mth discretized cell in the cathode, and air chan- 
nel, respectively. The entire solution procedure follows a space 
marching algorithm. At each axial position the transient system 
of equations are solved until a steady state solution is obtained. 
The initial condition at each axial position assumes the con- 
verged solution from the previous finite volume cell. Due to the 
elliptic nature of heat balance equation in porous electrodes and 
the conductive terms in the solid regions, an outer iteration loop 
is formed around the marching algorithm. The equation system 
is solved using the differential algebraic equation (DAE) solver 
LIMEX [23]. The entire solution converges in few passes. 


4. Results and discussion 
The effect of various parameters on cell performance is stud- 


ied systematically. The study covers the effect of (a) air flow 
rate, (b) anode thickness, (c) steam to carbon ratio of the fuel, 


(d) specific area available for surface reactions, and (e) extend of 
pre-reforming on cell efficiency and power density. In all cases, 
computations are performed for adiabatic as well as isother- 
mal operations. For adiabatic calculations the outer interconnect 
walls are assumed to be adiabatic, i.e., no heat transfer over the 
interconnects. All isothermal calculations are carried out for a 
constant temperature of 800 °C. Furthermore, in all cases the cell 
is assumed to operate at a constant voltage of 0.7 V. Table 1 lists 
conditions used for the various cases studied and the membrane 
electrode assembly (MEA); parameters are given in Table 2. 
In all cases, the fuel stream consisting of 40 vol.% CHaand 60 


Table 2 
MEA parameters for all simulations 
Parameters Values Units 
Anode 
Thickness (la) 0.75 mm 
Width (Wa) 1.00 mm 
Average pore radius (rp) 0.50 pm 
Average particle diameter (dp) 2.50 pm 
Specific area (As) 1025.00 cm7! 
Porosity (€) 0.35 
Tortuosity (Tt) 3.80 
Charge transfer coefficient (8a) 0.50 
Specific heat (C pr) 450.00 Jkg7! K7! 
Density 3310.00 kg m~? 
Thermal conductivity (k) 1.86 Jm! s7! K7! 
Electrolyte 
Thickness (le) 25.00 pm 
Width (W.) 1.00 mm 
Specific heat (Cpe) 470.00 Jkg-! K7! 
Density 5160.00 kgm? 
Thermal conductivity (k) 2.16 Im!s=!K=! 
Cathode 
Thickness (le) 30.00 pm 
Width (W.) 1.00 mm 
Average pore radius (rp) 0.50 pm 
Average particle diameter (dp) 2.50 pm 
Porosity (€) 0.35 
Tortuosity (7) 3.80 
Charge transfer coefficient (ß,) 0.50 
Specific heat (Cp) 430.00 Jkg-! K7! 
Density 3030.00 kg m~? 
Thermal conductivity (k) 5.84 Jm! s7! K7! 
Interconnect 
Thickness (4) 300.00 pm 
Specific heat (Cpr) 550.00 Jkg-! K7! 
Density 3030.00 kg m~? 
Thermal conductivity (k) 20.00 Jm! s7! K7! 
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Fig. 2. Temperature profile within flow channels and interconnects for a cell 
operating at 0.7 V. Air stream is assumed to enter the cathode channel at 650°C. 
The cell is assumed to operate under adiabatic conditions. 


vol.% H20 is assumed to enter the cell at 800° C with a velocity 
of 30cms-!. 

The temperature distributions in the fuel/air channels and 
the interconnects are shown in Fig. 2. Interestingly, the anode 
side interconnects are at lower temperature than the fuel temper- 
ature; nonetheless, the cathode side interconnects are at higher 
temperature than the air temperature. The increased temperature 
of cathode side interconnect compared to the air temperature is 
primarily due to the heat radiation from the cathode electrode 
to the interconnect. Near the inlet, the fuel stream transfers heat 
to the relatively cold air entering the cathode, further heat is 
consumed by the endothermic reforming reactions. 

The temperature distributions within the anode and cathode 
are presented in Fig. 3(a) and (b). The decrease of temperature 
near the fuel inlet due to reforming also affects fuel utilization. 
For direct internal reforming longer cells are required to achieve 
better fuel utilization. The fuel and air utilization for the case 
shown in the Fig. 3(a) and (b) are 57.5% and 24%, respectively. 
The temperature increase downstream the channel is primar- 
ily due to the exothermic cell reactions; overpotential losses, 
shift reactions, and the charge transfer reactions are exother- 
mic. There is a general notion that internal reforming leads to 
temperature drop near the cell inlet. However, the position of 
the reforming zone is determined by the inlet fuel velocity/flow 
rate. Fig. 4 demonstrates the influence of fuel inlet velocity on 
the position of reforming zone. In the case of 30 cms! the 
reforming zone is 2cm down the fuel inlet. However, for the 
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Fig.3. Temperature distribution within the anode and cathode for a cell operating 
at 0.7 V under adiabatic conditions. Air stream is assumed to enter the cathode 
channel at 650 °C: (a) anode and (b) cathode. 
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Fig. 4. Influence of the inlet velocity on the position of reforming zone. 


case of 10cms~! the reforming zone if 5mm down the fuel 
inlet. 

For an adiabatic case, the species distributions within the 
fuel channel and the anode are shown in Fig. 5. The top panel 
displays the species distribution along the fuel channel, while the 
drop down panels depict the profiles across the anode thickness 
at selected axial positions. The reforming chemistry within the 
cell produces Hy and CO by consuming CH4 and H20. Although 
the concentration of H20 initially decreases due to reforming, 
the electrochemical production of H20 leads to further increase 
downstream the cell. The concentration of Hz gradually drops 
as a result of electrochemical consumption. Although, CO can 
participate in charge transfer reactions, in the analysis carried 
out here, H2 is considered to be the only electrochemically active 
species. Therefore, any CO produced participates either in shift 
reaction or in Bouduard reaction. 

Fig. 6 displays the surface coverages of surface adsorbed 
species. The bottom panel shows the coverages along the three- 
phase boundary. The surface coverage of oxygen decreases 
initially and increases further down the length of the reactor, fol- 
lowing the trend of H20 concentration in the gas-phase. This is 
an indication that the surface coverage of oxygen results from the 
dissociative adsorption of H20. Most of the surface is covered 
by hydrogen and CO. The hydrogen coverage near the three- 
phase boundary initially increases and then drops. The increase 
in hydrogen coverage near the inlet is as a result of low elec- 
trochemical consumption due to the temperature drop near the 
inlet. This conclusion is supported by the profile of current den- 
sity along the length of the channel, where the current density 
decreases as a result of temperature decrease (Fig. 7). Though 
the mechanism is ineffective in predicting the carbon depo- 
sition quantitatively, qualitative agreement with experimental 
observations has been reported previously [7]. 

Though the results show negligible carbon deposition it 
should be kept in mind that the chemistry mechanism used 
here is not particularly validated for carbon deposition con- 
ditions. Downstream in the channel, the free Ni surface area 
increases possibly due to increasing desorption rate with increas- 
ing temperatures and due to the consumption of CO and H2 
by shift and electrochemical reactions. The drop up panels 
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Fig. 5. Species profiles within the fuel channel and anode. The drop down panels shows the species profiles across the anode thickness at various axial positions. 
The anode considered is 750 pm thick and the inlet fuel enters at 800°C, while air is assumed to enter at 650°C. 
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Fig. 6. Surface coverages of major adsorbed species. The drop down panels shows the species profiles across the anode thickness at various axial positions. The anode 
considered is 750 um thick and the inlet fuel enters at 800 °C, while air is assumed to enter at 650°C. The bottom panel shows the coverages along the three-phase 
boundary and the top panels display the coverages across the anode thickness. 
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Fig. 7. Current density as a function of the axial position for various anode 
thicknesses. Air enters the cathode channel at 650°C, X =1. 


present the coverages across the anode thickness at selected axial 
positions. 


4.1. Air flow rate 


Maintaining high air flow rates to keep the cell temperature 
within the limits of durable electrode performance is essen- 
tial for long term operation of SOFCs. However, the amount 
of excess air required depends on the inlet fuel composition. 
For example, an internally reforming cell requires less excess 
air than a cell operating on pure hydrogen fuel. However, it is 
well known that high operating temperature results in better cell 
performance, and hence an inordinate excess air flow rate can 
dramatically reduce the cell performance by lowering the aver- 
age cell temperature. Therefore, choosing the amount of excess 
air, resulting in an optimum balance between cell performance 
and cell durability is critical. 

In the calculations carried out here, excess air is defined in 
terms of air number A defined as 


A= BE. (30) 
Afuel Lmin 
where 
Omin 
Lin = : 31 
min x, (31) 


In the above equations, fuel is the molar flow rate of the 
fuel stream, Omin is the minimum O7 concentration required 
for the complete oxidation of the inlet fuel and [X 10, is the 
concentration of O> in the cathode stream. 

Fig. 8 shows the efficiency and power density as a function 
of air number (Eq. (30)). Fuel and air streams are assumed to 
enter the cell at 800 and 650 °C, respectively. It is quite apparent 
from Fig. 8 that both efficiency and power density decrease with 
increasing air flow rate, because in the case of adiabatic operation 
an increase in air flow rate decreases the overall cell temperature 
resulting in reduced average current density and hence decreased 
fuel utilization. The gradient as well as the maximum and the 
minimum temperature of the solid are shown in Fig. 9. 

However, the trend reverses when the cell is operated isother- 
mally (Fig. 10). In this case both efficiency and power density 
increase because of the reduced air depletion with an increasing 
air flow rate. Furthermore, isothermal operation results in higher 
efficiency and power density compared to non-isothermal case. 
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Fig. 8. Efficiency and power density for adiabatic condition as a function of air 
number. Cathode inlet stream (air) enters at 650°C. 
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Fig. 9. Maximum, minimum, and temperature gradient of the solid for varying 
air flow rate to fuel flow rate; air number according to Eq. (30). 


The maximum efficiency for non-isothermal operation is 44%, 
while isothermal operation results in a maximum efficiency of 
59%. It should be noticed that under constant operating potential 
(0.7 V) and inlet fuel compositions, the vital factor governing 
power density and efficiency is the fuel utilization defined by 
Eq. (26). Therefore, in an operating stack each cell can lead 
to different levels of fuel utilization and hence different effi- 
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Fig. 10. Efficiency and power density for isothermal condition (800°C) as a 
function of air number. 
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Fig. 11. Effect of anode thickness on efficiency and power density for a cell 
operating under adiabatic conditions. Cathode inlet is assumed to be air at 
650°C. 


ciency levels depending on the temperature distribution, which 
strongly depends on the cell’s position within the stack and the 
temperatures of the adjacent cells. 


4.2. Anode thickness 


Though SOFC can be either anode, electrolyte, or cathode 
supported, in the case of cells running on hydrocarbon fuels, 
anode supported cells may be preferable due to the potential for 
internal reforming. However, finding the optimal anode thick- 
ness required to support the cell mechanically and to achieve the 
desired level of internal reforming and optimal cell performance 
is a rather difficult task. 

Fig. 11 displays the influence of the anode thickness on effi- 
ciency and power density for a cell operating under adiabatic 
conditions. Quite surprisingly the model predicts an optimal 
anode thickness at ~0.5 mm maximizing the efficiency and 
power density. In every case, the fuel has to undergo internal 
reforming to produce H2 and CO which further participate in the 
charge transfer reactions at the three-phase boundary. However, 
in the case of thin anodes the short diffusion path (~0.2—0.4 mm) 
available for the fuel before reaching the three-phase boundary 
limits the amount of H2 and CO produced by internal reforming 
and, hence, leads to a low average current density, efficiency, and 
power density. As the thickness increases, the possibility for the 
fuel to undergo reforming to produce H2 and CO also increases. 
However, the higher extend of internal reforming achieved from 
the longer diffusion path also leads to larger temperature drops 
close to the channel inlet and hence results in decreased perfor- 
mance. Thus, there turns out to be an optimal anode thickness 
where the efficiency and power density reaches a maximum. 
However, it should be noticed that the optimum thickness can 
vary with the operating conditions, such as fuel composition, 
inlet temperature, velocity, catalyst loading, etc. 

Fig. 12 presents the efficiency and power density as a function 
of anode thickness for the case of a cell operated isothermally. 
In this case, both efficiency and power density increase with 
increasing anode thickness, however, reaches an asymptotic 
behavior at higher anode thicknesses. This clearly indicates that 
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Fig. 12. Effect of anode thickness on efficiency and power density for a cell 
operating under isothermal conditions. 


even in the case of isothermal operation increasing the anode 
thickness beyond the optimal value does not bring any fur- 
ther increase in performance. Nevertheless, isothermal operation 
results in better performance than adiabatic operation. A maxi- 
mum efficiency of 59% is achievable with isothermal operation, 
while the maximum possible efficiency in the case of adiabatic 
operation is 45%. 


4.3. Steam to carbon ratio 


It is quite well known that Ni based anodes are prone to cok- 
ing while running on hydrocarbon fuels [24—26]. Therefore, it is 
quite important to have enough steam present in the inlet fuel to 
suppress coking. Figs. 13 and 14 respectively present the influ- 
ence of steam to carbon (s/c) ratio on the performance of an 
adiabatically and isothermally operating cell. In both cases, the 
efficiency improves with increasing steam content in the inlet 
fuel stream, while the power density decreases with increasing 
steam content. The decrease in power density is due to the fuel 
dilution effect; diluted fuel results in a lower average current den- 
sity and for the cases studied here power density solely depends 
on the current density due to the constant operating voltage. 
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Fig. 13. Efficiency and power density for varying steam to carbon (s/c) ratio in 
the inlet fuel stream for a cell operating adiabatically. For all s/c ratio the fuel is 
assumed to enter at 800° C and air at 650°C. 
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Fig. 14. Efficiency and power density for varying steam to carbon ratio (s/c) in 
the inlet fuel stream for a cell operating isothermally at 800°C. For call cases, 
the cathode stream is assumed to be air. 


Though the average current density decreases, the efficiency 
increases with increasing steam content. This is essentially 
counter intuitive. One would expect decreasing efficiency with 
increasing steam dilution, which is certainly true for the case of 
maximum efficiency analysis. However, one should notice that 
the analysis carried out here is for a particular cell geometry, 
and the fuel utilization for a particular cell depends on the cell 
dimensions. From the definition of efficiency (Eq. (25)) it is evi- 
dent that for a given cell potential, efficiency strongly depends 
on the fuel utilization, because the thermoneutral voltage Ein is 
always ~1.2 V. At low s/c ratios the stoichiometric current is 
much higher than the average current delivered by the cell lead- 
ing to low fuel utilization and hence low efficiency. However, as 
the fuel dilution increases the stoichiometric current approaches 
the delivered current resulting in increased fuel utilization and 
hence efficiency. Fig. 15 shows the stoichiometric current and 
average current for the case of a cell operated at 800 °C. 

Here, like in the cases mentioned above isothermal operation 
results in improved cell performance. In the case of higher s/c 
ratio a higher efficiency resulting from the higher fuel utilization 
leads to reduced power density, which would require large cell 
areas to draw the required volumetric power density. On the 
other hand, operating conditions providing high power density 
result in a very low efficiency. Therefore, a judicious choice of 
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Fig. 15. Average current and stoichiometric current for varying steam to carbon 
ratio. 
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Fig. 16. Effect of specific area on efficiency and power density for cell operating 
adiabatically. Inlet fuel enters at 800 °C, and the air stream is assumed to enter 
the cathode at 650 °C. 


steam dilution has to be made to optimize power density and 
efficiency. 


4.4. Specific area 


Sufficient catalytically active surface area is needed for het- 
erogeneous reactions. In literature, the amount of catalyst used 
in the anode usually is expressed in terms of weight percentage 
(wt%). However, from chemisorption measurements of sample 
anodes, it is possible to retrieve information on the available 
active surface area. A large catalytically active surface area 
results from higher catalyst loading, better catalyst dispersion, 
and better microstructural properties of the electrode. Fig. 16 
shows the effect of specific area (area/volume) on the cell per- 
formance for a cell operated adiabatically. Both efficiency and 
power density are found to attain a maximum value at an avail- 
able specific area of ~1x 10° m7!. Large amount of catalysts 
can lead to excessive reforming, thereby decreasing the tem- 
perature near the inlet of the cell, resulting in a decreasing 
performance. A poor availability of catalytically active area 
results in lower production of synthesis gas participating in 
the charge transfer reactions, and hence resulting in a lower 
average current density. However, in the case of isothermal oper- 
ation, both power density and efficiency are found to increase 
with increasing specific area (Fig. 17) and reaches an asymp- 
totic behavior at higher values. Therefore, even in the case of 
isothermal operation, there is an optimal value for specific area 
beyond which there is no additional benefit by increasing the 
catalyst loading or dispersion. Nevertheless, isothermal opera- 
tion results in better efficiency and power density compared to 
adiabatic operation. 


4.5. Pre-reforming 


Though the focus of this paper is on direct internal reform- 
ing, the existing applications use some extend of pre-reformed 
fuel. Therefore, a systematic study to understand the influence of 
non-reformed and pre-reformed fuels on cell efficiency is car- 
ried out. In the case considered here the non-reformed fuel is 
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Fig. 17. Effect of specific area on efficiency and power density under isothermal 
condition (800 °C). 


assumed to consist of 60 vol.% CH4 and 40 vol.% H20; 40% 
pre-reforming means 40% CH4 converted in the pre-reformer. It 
is expected that direct internal reforming can result in reduced 
cost and increased overall efficiency of the system and can also 
help to resolve the temperature distribution problem in the cell. 
However, it is quite convincing from Fig. 18 that the efficiency 
of the fuel cell is higher for pre-reformed fuel. Both efficiency 
and power density increase with extent of pre-reforming for both 
adiabatic and isothermal case (Fig. 19). Unlike other cases, adi- 
abatic operation results in better performance of the cell with 
increasing extent of pre-reforming. For example, with 60% pre- 
reforming, adiabatic operation results in ~45% efficiency, while 
isothermal operation results in ~42% efficiency. However, with 
the lowest extend of pre-reforming considered here (20%), adia- 
batic operation results in lower efficiency and power density than 
isothermal operation. With increasing extent of pre-reforming, 
the H> content of inlet fuel increases, resulting in higher cell tem- 
peratures leading to higher performance, while the non-reformed 
fuel results in internal reforming leading to temperature drop 
and hence lower power density and efficiency. Furthermore, due 
to the short reactor length considered here, the fuel utilization 
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Fig. 18. Effect of pre-reforming on efficiency and power density under adiabatic 
condition. In call cases the pre-reformed fuel is assumed to enter at 800° C and 
air at 600 °C. The non-reformed fuel is assumed to be 60 vol.% CHy and 40 vol.% 
H20. 
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Fig. 19. Effect of pre-reforming on efficiency and power density under isother- 
mal condition (800 °C). The non-reformed fuel is assumed to be 60 vol.% CH4 
and 40 vol.% H20. 


efficiency turns out to be lower for non-reformed fuel. This also 
leads to a lower overall efficiency. However, it should be noticed 
that with increasing levels of pre-reforming, increasing amounts 
of excess air to keep the cell temperature within the limits would 
be required, and the system should be designed to utilize the 
waste heat produced by the cell. 


5. Conclusions 


Based on detailed models of transport and chemistry we have 
systematically analyzed the effect of various parameters on cell 
performance for a particular SOFC geometry. The vital factor 
for optimal cell performance has been found to be the resulting 
temperature profile within the cell. Judicious choice of steam 
dilution of the inlet fuel has to be made, that can well suppress the 
coking, while achieving reasonable efficiency and power den- 
sity. As aremarkable outcome, excessive amount of catalysts and 
highly thick anodes have been found to have detrimental effect 
on cell performance in the case of direct internal reforming under 
adiabatic conditions. Even for the case of isothermal operation 
optimal values exists for anode thickness and active specific 
area. Furthermore, though direct internal reforming is known 
to increase the overall efficiency of the system, the efficiency 
of the fuel cell itself is found to be higher for the case of pre- 
reformed fuel. In all cases, except for the case of pre-reformed 
fuel, isothermal conditions resulted in better cell performance. 
Therefore, in the case of direct internal reforming, it is quite 
important to have the cell designed in a way to achieve more or 
less isothermal operation of the stack, which can result in higher 
performance. 
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